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Large spin-phonon coupling and magnetically-induced phonon anisotropy in SrMC>3 

perovskites (M=V,Cr,Mn,Fe,Co) 

Jun Hee Leqj and Karin M. Rabe 

Department of Physics and Astronomy, Rutgers University, Piscataway, New Jersey 08854-8019, USA 

First-principles calculations reveal large zone-center spin-phonon coupling and magnetically- 
driven phonon anisotropy in cubic perovskites SrMC>3 (M=V,Cr,Mn,Fe,Co). In particular, the 
frequency and splitting of the polar Slater mode is found to depend strongly on magnetic ordering. 
The coupling is parameterized in a crystal-structure-dependent Heisenberg model and its main fea- 
tures seen to arise from the Goodenough-Kanamori rules. This coupling can be expected to produce 
distinct low-energy alternative phases, resulting in a rich variety of coupled magnetic, structural 
and electronic phase transitions driven by temperature, stress, electric field and cation substitution. 

PACS numbers: 75.80.+q, 63.20.-c, 75.10.Hk, 77.80.-c 
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Multiferroic phases, which exhibit multiple order- 
ings such as spontaneous polarization, magnetization, 
and strain, are currently the subject of much active 
research^. The main interest lies not in the coexistence 
of the orderings in the same material but rather, in the 
couplings between the orderings. In combination with 
direct control of the orderings by macroscopic fields and 
stresses, these couplings lead to physical properties in- 
cluding linear magnetoelectric response, magnetodielec- 
tric response and cross-switching of magnetization by 
electric fields or polarization by magnetic fields. Fur- 
thermore, multiferroic systems can exhibit coupled phase 
transitions in which the system transforms from the bulk 
phase to an alternative low-energy phase with distinct or- 
dering of multiple types, with functional behavior at the 
phase boundary arising from switching between the two 
phases. This behavior offers both fundamental scientific 
interest and the possibility of novel multifunctionality for 
technological applications such as sensors, actuators and 
information storage. 

One major challenge is that magnetostructural cou- 
pling is generally quite small. Magnetocrystallinc 
anisotropy, which measures the coupling of the direction 
of the spin magnetization to the lattice, depends on the 
spin-orbit interaction and is extremely small for all but 
the highest atomic number elements. Spontaneous po- 
larization is invariant under time reversal symmetry, the 
principal symmetry broken by spontaneous magnetiza- 
tion, and conversely, spins are invariant under inversion, 
the principal symmetry broken by spontaneous polariza- 
tion. In improper multiferroics such as TT^MnOg^, com- 
plex magnetic ordering breaks inversion symmetry and 
produces polarization, but the induced polarization is 
tiny. Moreover, the frustration needed to produce this 
magnetic ordering leads to low transition temperatures, 
while higher operating temperatures are needed for prac- 
tical applications. 

On the other hand, certain changes in interatomic dis- 
tances and bond angles can have a large effect on ex- 
change couplings, leading to changes in magnetic order- 
ing energies. Systems exhibiting this magnetostructural 
coupling can be efficiently identified by first-principles 



calculations^— of the change in phonon frequency with 
magnetic order, referred to as spin-phonon coupling. Sys- 
tems where the energy of these structural changes is com- 
parable to the change in magnetic energy can have large 
responses to applied fields and the possibility of cou- 
pled phase transitions. For example, in EuTiOg^ and 
SrMnOyi, it was found that as a result of spin-phonon 
coupling, a ferromagnetic phase with a ferroelectric polar 
distortion is competitive in energy with the bulk antifcr- 
romagnetic paraelectric phase, leading to a coupled phase 
transition driven by epitaxial strain. 

In this paper, we report a systematic first- 
principles study of the spin-phonon coupling of zone- 
center polar phonons in the cubic perovskites SrA/03 
(M=V,Cr,Mn,Fe,Co), with results for non-magnetic 
SrTiOa included for comparison. For all magnetic ma- 
terials studied, the Slater-type phonon mode is found 
to exhibit the largest spin-phonon coupling, as mea- 
sured by the frequency shift between the FM and G- 
AFM ordered phases. Further study of SrMn03 and 
SrCo03, which show the largest shifts, reveals pro- 
nounced magnetically-induced phonon anisotropy. This 
large spin-phonon coupling is parameterized with a 
crystal-structure-depcndent Heisenberg model and ana- 
lyzed with the help of the Goodenough-Kanamori rules^. 
Experimentally observable consequences of the spin- 
phonon coupling are discussed, in particular the possibil- 
ity of low-energy alternative phases and coupled phase 
transitions as a function of applied fields, stress, and 
cation substitution. 

First-principles calculations were performed using 
density-functional theory within the generalized gradi- 
ent approximation GGA+U method^ with the Pcrdcw- 
Becke-Erzenhof parameterization^ as implemented in 
the Vienna Ab Initio Simulation Package (VASP- 
4.6) 12 ' 13 . We use the Dudarev^ implementation with on- 
site Coulomb interaction £7=2.5 eV and on-site exchange 
interaction Jjj=1.0 eV to treat the localized 3d electron 
states. The value of U for Mn was adjusted to 2.7 eV to 
fit the experimental magnetic moment of 2.6±0.2 ug for 
SrMnOgiii. For Ti, we use [7 = 0. The projector aug- 
mented wave (PAW) potentials^ explicitly include 10 va- 



TABLE I: Cubic perovskite lattice constant a^ (in A), mag- 
netic ordering type and band gap (in eV), compared with 
available experimental data (in parentheses). The lattice con- 
stant was calculated for the computed magnetic type given in 
the table, except for that of SrCr03, which was computed for 
G-type ordering and compared with measurement in the para- 
magnetic phase. If the system is metallic, the gap is reported 
as m. The relative energies AE of alternative magnetic order- 
ings specified by the letter in parentheses (in meV/f.u.) are 
computed at the given computed lattice constant. The ex- 
perimentally determined ordering temperatures T c (in K) are 
included. Experimental data for SrVC>3 is from Refs. Jlall9l . 
for SrCrOs from Ref. J2I for SrMn0 3 from RefsJpEJT for 
SrFeOs from Refs. JMH for SrCoOs from Ref. Ijand for 
SrTiOs from Refs. [MM- 
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gap 



AE 



T c 



SrV0 3 


3.88(3.85) 


F(SDW) 


m(m) 


37(G) 


85 K 


SrCrOs 


3.86(3.82) 


C(C) 


m(m) 


7(G),53(F) 


40 K 


SrMn0 3 


3.85(3.81) 


G(G) 


0.45(i) 


86(F) 


233 K 


SrFe0 3 


3.89(3.85) 


F(hel) 


m(m) 


220(G) 


134 K 


SrCoOg 


3.84(3.83) 


F(F) 


m 


258(G) 


305 K 



SrTiOs 3.95(3.91) N(N) 1.8(3.3) 



lence electrons for Sr (4s 2 4p 6 5s 2 ), 6 for oxygen (2s 2 2p 4 ), 
13 for V (3s 2 3p 6 3d 4 4s 1 ), 12 for Cr (3p 6 3d 5 4s 1 ), 13 for Mn 
(3p 6 3dHs 2 ), 14 for Fe fip^dHs 1 ), and 9 for Co (ScPis 1 ). 
The zone-center phonon frequencies of the ideal cubic 
perovskite reference structures with various magnetic or- 
dcrings were computed using the frozen phonon method 
in supercells with y/2 x y/2 x \[2, \[2 x \/2 x 1, 1x1x2, 
and lxlxl primitive perovskite cells for G-AFM, G- 
AFM, A-AFM, and FM, respectively. Spin-orbit cou- 
pling was not included, so that the phonon wavevectors 
for the primitive perovskite cell remain good quantum 
numbers. We obtained the modes at q = by uniformly 
displacing all atoms in related by translation symmetries 
of the primitive perovskite structure. For each super- 
cell, Monkhorst-Pack (M-P) /c-point meshes 1 ^ were cho- 
sen to obtain phonon frequencies converged to within a 
few cm -1 . 

The computed ground state cubic lattice parameters 
and magnetic orderings are reported in Table U and com- 
pared with experimental information. The lowest energy 
magnetic ordering is correctly computed in all cases ex- 
cept SrVOj 1 ^ and SrFeO^, which are observed to have 
ferromagnetic spin density wave and complex helical or- 
dering, respectively, not included in the first-principles 
analysis. The computed relative energies of higher- 
energy magnetic orderings are reported and are seen 
to correlate roughly with the experimentally-observed 
magnetic-ordering temperatures given in Table HI The 
computed lattice constants are larger than experimen- 
tal values by about 1%, as typical for GGA calculations 
of oxides. Consistent with experiment, all the magnetic 
compounds are found to be metallic except for SrMn03, 
which has a computed gap of 0.45 eV, comparable to that 



TABLE II: Computed IR-active polar phonon frequencies 
(cm -1 ) of G-AFM and FM states in cubic SrM0 3 perovskites. 
The calculations were done with computed lattice constants 
given in Tabled] 



SrTiOs {(f) 
non-mag. 



SrVOs (d 1 ) 
G-AF FM 



SrCrOs (d 2 ) 
G-AF FM 



Slater 


130 i 


329 


290 


287 


187 


Last 


144 


148 


147 


156 


145 


Axe 


508 


541 


546 


525 


512 



SrMnOa (d 3 ) 
G-AF FM 



SrFeOs (d 4 ) 
G-AF FM 



SrCoOa (d 5 ) 
G-AF FM 



Slater 120 122 i 
Last 165 165 

Axe 475 486 



248 231 
148 150 
506 505 



176 i 203 
155 154 
508 498 



found in a previous study 2 ^. 

In Table [TT] we report the computed zone-center polar 
phonon frequencies for FM and G-AFM ordering. The 
phonons are labeled according to the character of the 
eigenmode. The Slater mode 2 ^ involves the oscillation 
of the B-cation against the oxygen-octahedron network. 
This is the soft mode which is closely associated with 
the ferroelectric instability in perovskites such as BaTiOs 
and KNb03. The Last mode 2 ^ involves the oscillation of 
the A-cation against the oxygen-octahedron network. It 
is this mode which is the soft mode for A-site ferroelec- 
tricity in perovskites such as PbTiC>3 or BiFe03. The 
Axe mode^S, at the highest frequency, corresponds to 
distortion of oxygen octahedra. 

While the Last and the Axe mode frequencies are, as 
would be expected from their eigenmode character, rel- 
atively insensitive to the -B-site <i-occupancy and mag- 
netic ordering, the Slater mode shows a dramatic shift 
with change in d-occupation from compound to com- 
pound, and with change in magnetic ordering for a given 
compound. The Slater modes are highest in frequency 
for M=V (d 1 ), M=Cr (d 2 ) and M=Fe (d 4 ); with cubic 
symmetry ti g -e g splitting, these cases involve partially 
occupied sublevels. For M—Mn (d 3 ) and M=Co (d 5 ), 
the Slater modes are substantially lower in frequency. 
The frequency shifts of the Slater mode between FM and 
G-AFM ordering are more than 100 enr 1 for SrCrOs, 
SrMn0 3 and SrCo0 3 . In SrV0 3 and SrFe0 3 , though the 
shift is smaller than for the other compounds, it is still 
significantly larger than for the other two modes. The 
low frequencies and the large magnetic ordering shifts in 
SrMnC>3 and SrCoC>3 lead, in fact, to ferroelectric insta- 
bilities in the FM and G-AFM orderings, respectively. 
However, in each case these are not the lowest energy 
magnetic orderings, and thus the ground state structure 
does not exhibit a polar distortion. Next, we computed 
the polar phonon frequencies with two lower-symmetry 
magnetic orderings for SrMn03 and SrCo03, the two sys- 
tems with the largest frequency shift with magnetic or- 
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FIG. 1: (Color online) Computed frequencies for the Slater 
mode(s) for various magnetic orderings of cubic SrMnOs, with 
results for SrCo03 given in parentheses. Arrows represent 
atomic displacements. Black and shaded circles represent up 
and down spins on the magnetic _B-site in ABO3 perovskites; 
small open circles represent oxygen atoms. Frequencies are 
given in cm - . 



dering. Specifically, we considered C-type, with AFM 
ordering in the xy plane and FM ordering along the nor- 
mal direction, giving periodicity q = (it /a, 7r/a, 0), and 
A-type, with FM ordered planes ordered antifcrromag- 
netically along the normal direction, giving periodicity 
q= (0,0, n/a). 

The magnetic ordering and the frequencies of the as- 
sociated Slater modes are shown in Figure [T] The sym- 
metry lowering due to the C- and A-AFM orderings lifts 
the degeneracy of the Slater mode in each case by a con- 
siderable amount. This frequency splitting, which occurs 
even without any symmetry-lowering relaxation of the 
crystal structure, is referred to as magnetically-induced 
phonon anisotrop y 31 ' 32 . In C(xy)-AFM SrMnOs, the 
Slater mode with polarization along z is much stiffer than 
that with polarization in the plane, which is unstable; the 
splitting is reversed in SrCoOa, with both modes unsta- 
ble. In contrast, in A(z)-AFM SrMnOs both modes are 
unstable though the ferroelectric instability along the z- 



direction is favored, while in SrCoOa the Slater mode 
with polarization along z is much stiffer than that with 
polarization in the plane, which has a very low frequency. 

Each magnetic ordering and atomic displacement pat- 
tern can be characterized by the B spins and atomic dis- 
placements of B and oxygen relative to the line B-O-B: 
same spin and transverse displacement (FT), opposite 
spin and transverse displacement (AT), same spin and 
longitudinal displacement (FL) and opposite spin and 
longitudinal displacement (AL). For SrMnOs, low fre- 
quency is primarily associated with the FT bonds, and 
to a lesser extent the AL bonds, which predominate in 
A(z) magnetic ordering with displacements along the z 
direction. On the other hand, for SrCoOs low frequency 
is associated with the AT bonds. 

For a precise quantitative description of the spin- 
phonon coupling, we use the model parameterization of 
Ref. [3ll The magnetic interactions are described by a 
Heisenberg model for the spins on the B sit o 33 ' 34 , with 
exchange couplings determined by the computed energy 
differences between selected magnetic orderings. The 
coupling to the lattice originates from the dependence 
of the exchange coupling parameters on the atomic posi- 
tions. Thus, changes in atomic displacement patterns re- 
sult in changes in magnetic ordering energies; conversely, 
changes in magnetic orderings result in changes in the en- 
ergetics of atomic displacements, in particular the inter- 
atomic force constants, producing the spin-phonon cou- 
pling. 

For the cubic magnetic perovskite compounds con- 
sidered here, the model energy is written E = Eq + 
E p h + E S pi n . Eq is the energy of the paramagnetic 
cubic perovskite structure. The phonon term E p h = 
5 Erar'a' K TaT i a 'U Ta u T > a i is written here to include only 
primitive-perovskite cell zone-center modes, so that t in- 
dexes the five atoms in the unit cell and a indicates 
the Cartesian direction of the sublattice displacement u T . 
= — J2i^j JijSi ' Sj (5=3/2), where i and j index 



E 



spin 

the B site spins, is the Heisenberg model for the magnetic 
ordering energy. The exchange couplings arc included up 
to third neighbors to reproduce the relative magnetic or- 
dering energies in the cubic reference structure u T ~0. 
For SrMn0 3 , this yields Ji=-3.1 mcV/f.u., J 2 =-0.16 
meV/f.u. and J 3 =0.0027 meV/f.u. For SrCo0 3 Ji=8.8 
meV/f.u., J 2 =0.59 meV/f.u. and J 3 =0.69 meV/f.u. The 
decay with interspin spacing is considerably less rapid 
in SrCo03 than in SrMnOs, as expected for a metallic 
system. In the following, we will consider the atomic- 
displacement dependence only for the first two neighbor 
interactions. 

To investigate the spin-phonon coupling, we expand 
the model energy in the sublattice displacements u T . For 
spin configurations of sufficiently high symmetry, includ- 
ing the five configurations considered here, the first order 
terms vanish. At second order, the full force constant 



TABLE III: J(' a (upper) and J'{ z (lower) of SrMn0 3 (left) and SrCo0 3 (right). Unit is in 1(T 4 eV/A 2 



T" 


Sr 


Mn 


0* 


Oy 


o z 


J" 


Sr 


Co 


o. 


O y 


o z 


Sr 


-61 


-14 


17 


70 


-9 


Sr 


-7 


-40 


66 


-17 


-5 


Mn 


-14 


2400 


-2000 


310 


-650 


Co 


-40 


-3300 


2400 


-170 


1100 


0, 


17 


-2000 


1500 


-43 


530 


0,- 


66 


2400 


-1700 


49 


-820 


°u 


70 


310 


-43 


-270 


-80 


o„ 


-17 


-170 


49 


83 


50 


o 2 


-9 


-650 


530 


-80 


210 


o z 


-5 


1100 


-820 


50 


-330 


J" 


Sr 


Mn 


o, 


o y 


0* 


T" 


Sr 


Co 


0* 


O y 


0, 


Sr 


-30 


-75 


26 


26 


54 


Sr 


27 


46 


-11 


-11 


-53 


Mn 


-7.") 


-1500 


220 


220 


1200 


Co 


46 


-1100 


250 


250 


540 


0,, 


26 


220 


-150 


-8 


-80 


o x 


-11 


250 


150 


19 


-410 


o y 


26 


220 


-8 


-150 


-80 


o„ 


-11 


250 


19 


150 


-410 


0, 


54 


1200 


-80 


-80 


-1100 


o, 


-53 


540 


-410 


-410 


330 



matrix K = d 2 E / du Ta du T i a i is given by 



K TC 



K T 



d 2 J tj 



*> Ta LL T ' q.' 



< Si ■ Si > 



By symmetry, the 15 x 15 



with </■'■_„ _/„, = 
matrices </"■ are block diagonal, with nonzero elements 
only for sublattice displacements u T and u T i along the 
same Cartesian direction. Furthermore, cubic symmetry 
relates the 5x5 block for displacements for one choice of 
ij and Cartesian displacement direction to the blocks for 
other choices of ij and Cartesian displacement direction, 
so that there are only four independent blocks: two for 
first neighbor exchange couplings (z-displacemcnt blocks 
J" x and J" z , with J" y related to J" x ,) and for next nearest 
neighbors (J^',, for displacements along z and exchange 
couplings within the xy plane, and J^± for displacements 
along z and exchange couplings to spins at greater or 
lower z). For the spin configurations considered, the force 
constant matrices have the same block diagonal form as 
the J" matrices and can be compactly expressed in terms 
of J" z-displaccmcnt blocks as follows: 



f K F (z) 1 




k G {z) 




K c (z) 




k c {x) 




K A (z) 




I K A (y) J 





-9/2 


-9/2 


-9/2 


-9 


+9/2 


+9/2 


+9/2 


-9 


+9/2 


+9/2 


-9/2 


-9 


+9/2 


-9/2 


+9/2 


+9 


-9/2 


-9/2 


+9/2 


-9 


-9/2 


+9/2 


-9/2 


+9 



-181 




r K i 


■18 




J lx 


■18 




J" 







T" 


■18 




J" 

2 II 







T" 



where K F , K G , K c( - x y\ K A ^ arc 5x5 blocks of the full 
force constant matrices. Using first-principles computa- 
tions for these matrices, we can determine the J" blocks 
and K, the force constant matrix for paramagnetic or- 
dering < Si ■ Sj >=0. 

Table Hill shows J" x and J[' z 5x5 z-displacement blocks 
for SrMnC>3 and SrCoOa. O^ and O y denote the equa- 
torial oxygens in the line B-O-B along x and the equa- 
torial oxygen in the line B-O-B along y, respectively. 



From these matrices, we find that the large spin-phonon 
coupling noted for SrMnOa and SrCoOa can be mainly 
attributed to the nearest neighbor exchange interaction 
submatrix J" x , specifically the transverse displacement 
of the B cations and oxygen atoms in the B-O-B bond 
along x. For a more detailed discussion, we diagonal- 
ized the J' lx matrices. Eigenvalues of the J" x block for 
SrMn0 3 are (0.42,-0.04,-0.01,0.005,0.000) with eigenvec- 
tor (0.00,0.76,-0.62,0.06,-0.21) for the largest eigenvalue. 
A phonon with this displacement pattern thus would ex- 
hibit the maximum spin-phonon coupling, explaining the 
large effect for the Slater mode. The fact that the four of 
the five eigenvalues are close to zero allows us to approxi- 
mate the atomic-displacemcnt-dcpcndcnt part of J\ x us- 
ing this single eigenvector: 

1 2 

AJ lx - 0.066 [(ll M n - UoJ + iKn ~ Woj] ■ 

The bond angle reduction represented by the first term 
makes the main contribution to the spin-phonon cou- 
pling, with the bond length change represented by the 
second term being a secondary contribution. 

For SrCoOs, the analysis is very similar. Eigen- 
values of the J[' x of SrCo0 3 are (-0.54,0.013,0.005,- 
0.003,0.000) with eigenvector (0.01,0.78,-0.56,0.03,-0.26) 
for the largest eigenvalue. Again, the other eigenvalues 
are close to zero, allowing the approximation 

1 2 

AJi^ - -0.093 [(uco - "oj + o ("Co - uoj] ■ 

As for SrMn03, the two terms represent the bond- angle 
reduction and bond-length change contributions, which 
for SrCoOs are negative rather than positive. 

The importance of the two terms in SrMn03 can 
be qualitatively understood using the Goodenough- 
Kanamori rules^. The magnetic exchange coupling J is 
given by J ~ £ 2 /(A — Jh) where t is the hopping integral, 
which is proportional to the effective wavefunction over- 
lap between the two cations through their shared oxygen, 



A is the energy difference between oxygen and B-cation 
orbitals and Jh is the Hund coupling. Reduction of the 
bond angle ((A6>) 2 ~ (ub — uoj 2 ) decreases the hopping 
integral t, which reduces the magnitude of the exchange 
coupling J. Also, displacement of the apical oxygens O z 
towards the B cation should shift the i3-cation and 2 
orbital energy level in opposite directions due to elec- 
trostatics; the resulting decrease in the magnitude of J 
arising from the change in A. For SrCoOs, ferromag- 
netism can be explained by analogy to the Zener double 
exchange (DE) mechanis m 35 ! 36 with high-spin Co d 6 ions 
and ligand holes of 1/3 per oxygen antiferromagnetically 
coupled to the Co spin, with the mobile holes favoring fer- 
romagnetic alignment of the cobalt spins 3 -^. If the d — p 
matrix element is decreased by a reduction in bond angle, 
or the change-transfer gap is increased by displacement 
of the apical oxygen, this will decrease the magnitude of 
the magnetic exchange coupling 3 ^. 

While, as discussed above, the Slater mode maximizes 
the spin-phonon coupling, the other polar modes show 
very little change with magnetic ordering] 6 ! 39 ! 40 The Axe 
mode consists mainly of 2 displacements, which reduce 
the B-O z bond length without changing the B-O x -B 
bond angles. The first term does not contribute and the 
change in J is much smaller, corresponding to a change 
only of A. The Last mode, which is a displacement of 
A-site cation with respect to the B-Oq network, does not 
change the relevant bond angles or bond lengths, and 
both terms are negligible. 

The large spin phonon coupling for the Slater mode 
found from first principles can be experimentally ob- 
served in bulk samples by changing the magnetic order- 
ing. In the paramagnetic phase, application of a mag- 
netic field slightly increases alignment of the spins, lead- 
ing to a phonon frequency shift that is rather small. For 
example, a mean field theory calculation for SrMn03 
gives a shift of approximately 0.1 cm" 1 with a field of 
5 T at T]sj. A much more dramatic shift can be observed 
by changing the temperature through the magnetic phase 
transition. As the magnetic ordering changes from para- 
magnetic above the magnetic ordering temperature to 
the low-temperature ordered phase, the change in the 
spin correlations produces a substantial shift in phonon 
frequency that is roughly half the G-AFM-FM splitting. 
Indeed, motivated by these first principles results, far IR 
measurements have been made on SrMnOg^ -, showing 
the expected hardening of the Slater mode on cooling 
through T N =250 K. 

Another readily measurable physical property is the 
magnetodielectric effect. For example, in La2MnNiOg, 
an applied magnetic field shifts the critical temperature 
of a phase transition at which the dielectric constant 
jumps by about 10%, producing a large magnetodielec- 



tric response between the two temperatures 4^ In addi- 
tion, magnetically-induced anisotropy could be directly 
observed by measuring phonon splittings in a material 
that has a C-AFM bulk phase, for example SrCrOs. In 
this material, a large coupling between the magnetic mo- 
ment and a breathing phonon in a hypothetical G-AFM 
phase has been noted from first principles^. 

In thin films and superlattices, the large spin-phonon 
coupling can produce a coupled magnetic-ferroelectric 
phase transition with epitaxial strain as polarization- 
strain coupling increases the instability of the polar mode 
for non-ground-state magnetic ordering (it should be 
noted that other perturbations, such as compositional 
substitution, could also be effective) . An epitaxial-strain- 
induced multifcrroic (ferromagnetic-ferroelectric) phase 
based on this mechanism was identified from first prin- 
ciples investigation in EuTiO^ and subsequently exper- 
imentally confirmed^ 3 -. First-principles investigation of 
the epitaxial-strain phase diagram of SrMn03 has shown 
an analogous multiferroic phase, with higher ordering 
temperature 6 -. In SrCo03, with a ferromagnetic metal- 
lic bulk state, the polarization-strain coupling increases 
the ferroelectric instability of the insulating antiferro- 
magnetic phase, driving a coupled magnetic-ferroelectric 
metal-insulator transition with epitaxial strain, which 
has been the subject of a first-principles study^i. These 
systems are expected to exhibit desirable functional be- 
havior at the phase boundary associated with switching 
between the two phases by applied fields or stresses. 

In summary, we have presented a first-principles inves- 
tigation of the spin-phonon coupling in the Sr ^4-site per- 
ovskite compounds SrA/0 3 (iU=V,Cr,Mn,Fc,Co). The 
coupling is largest for the Slater phonons, especially 
for SrMn03 and SrCo03 which show large frequency 
shifts with changes in magnetic ordering and magnet- 
ically induced phonon anisotropy. We show that this 
behavior can be qualitatively understood on the basis 
of Goodcnough-Kanamori rules. The computed spin- 
phonon coupling can be experimentally investigated di- 
rectly through measurement of the phonon frequency 
shift through the magnetic transition. Other conse- 
quences include the possibility of coupled phase tran- 
sitions with epitaxial strain or compositional substitu- 
tion and desirable functional behavior for systems at 
the phase boundary. Thus, computation of spin-phonon 
coupling can be a valuable screening tool in the first- 
principles design of novel functional and multifunctional 
materials. 
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